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ABSTRACT 


By  analyzing  CONREP  data  for  oiler  (AO)  and  ammunition  (AE)  ships 
replenishing  destroyers  (DD)  and  attack  carriers  (CVA)  with  regard  to 
day  and  night  operation  and  load  requirement,  it  is  shown  through  Chi- 
Square  goodness-of -f it  tests  that  alongside  time  data  can  be  fitted  to 
an  Erlang  or  exponential  distribution.   In  addition,  several  methods 
for  estimating  the  shape,  scale,  and  shift  parameters  for  gamma  and 
Erlang  distributions  are  presented. 

Employing  well  known  but  analytically  little  used  features  of  the 
replenishment  operation,  a  computer  simulation  model  is  then  formulated 
and  programmed  based  on  these  distributions,  and  thus  sensitive  to 
changes  in  load  requirements.   Using  the  simulation  as  an  experimental 
device,  an  example  is  run  to  demonstrate  its  use  in  estimating  the  time 
to  complete  the  operation  and  conducting  sensitivity  analyses  on  the 
various  load  requirements. 
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I.   INTRODUCTION 

The  primary  reason  the  United  States  Navy  is  capable  of  performing 
its  regularly  assigned  mission  as  well  as  assist  in  global  contingencies 
effecting  United  States'  interests,  is  the  ability  to  swiftly  resupply 
combatant  ships  almost  anywhere  in  the  world.   This  is  possible  through 
replenishment  at  sea  operations  which  render  ships  independent  of 
costly,  vulnerable  land  based  supply  facilities  while  effectively 
increasing  combatant  force  strength,  since  ships  are  able  to  remain 
on  or  near  their  assigned  stations  during  replenishment. 

Although  very  important,  replenishment  at  sea  is  a  hazardous 
operation.   During  replenishment,  participant  ships  are  extremely 
close  increasing  their  possibility  of  collision  and  vulnerability  to 
enemy  attack.   Even  though  vertical  replenishment  (VERTREP)  using 
helicopters  has  helped  to  reduce  the  risk  in  some  cases,  all  fuel 
transfers  such  as  Navy  Special  Fuel  Oil  (NFSO) ,  Aviation  Gasoline 
(AVGAS),  and  JP-5,  and  most  transfers  of  dry  material  such  as  ammuni- 
tion and  stores,  are  necessarily  done  by  the  more  dangerous  connected 
replenishment  (CONREP)  method.   Recognizing  the  importance  of  CONREP, 
the  U.S.  Navy  is  using  a  considerable  amount  of  its  resources  to 
increase  the  efficiency  and  reduce  the  danger  of  this  operation. 
Increased  training  of  personnel,  redesign  or  replenishment  ships  and 
equipment,  and  studies  utilizing  analytical  and  computer  simulation 
models  are  efforts  that  have  been  and  are  being  made  in  this  endeavor. 

Operations  Research  studies  of  the  replenishment  at  sea  operation 
have  been  very  few.   However,  McCoullough  l4 j  formulated  two  models  of 
a  simplified  underway  replenishment  operation.   One  an  analytical 
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infinite  cycle  model  using  multi-stage  queueing  theory  and  the  other 
a  non-cyclic  computer  simulation.   In  both  models  the  distribution  of 
service  time  at  each  replenishment  ship  was  assumed  to  be  exponential. 
Waggoner  [7~]   and  to  some  extent  Patterson  \j>]  also  used  the  assumption 
of  exponential  service  (i.e.  alongside)  times  in  replenishment  at  sea 
operations.   Besecker  \_l\   questioned  the  validity  of  this  assumption 
and  analyzing  CONREP  data,  showed  that  the  distribution  of  service  times 
can  be  fitted  with  a  gamma  distribution.   Categorizing  the  data  as  to 
participants  and  day  and  night  operation,  he  contended  that  the  data 
would  also  fit  the  more  useful  Erlang  distribution  in  each  case.   This 
paper  continues  the  investigation  of  replenishment  at  sea  operations 
along  the  lines  of  reference  \_1  \  . 
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II.   OBJECTIVES 

The  most  immediate  gain  in  CONREP  efficiency  can  be  experienced 
through  proper  planning  of  the  operation.   However,  studies  of  the 
operation  have  been  largely  ignored  in  the  planning  stage,  since  they 
lack  the  realism  necessary  to  gain  the  confidence  of  operational 
commanders.   Previous  models  (such  as  \W]      and  \i7~\  )  have  not  taken 
full  advantage  of  the  many  known  features  of  CONREP,  for  instance, 
arrivals  at  supply  ships  are  not  random  but  may  be  promulgated  in 
advance  by  the  Task  Group  Commander.   Also,  the  quantity  required  by 
each  combatant  ship  is  known  prior  to  the  operation  to  some  degree  of 
accuracy  rather  than  being  a  random  variable.   If  the  transit  time 
between  various  supply  ships  is  considered  deterministic  and  known, 
the  uncertainty  in  the  system  is  a  result  only  of  the  unknown  alongside 
times,  which  should  be  dependent  on  the  participants  and  the  loads 
transferred. 

The  first  objective  of  this  paper  is  to  develop  a  more  accurate 
estimation  of  alongside  times  for  various  combinations  of  CONREP 
participants  by  considering  a  given  load  requirement.   Such  conditional 
distributions  would  be  applicable  in  analytical  models,  in  computer 
simulations,  and  possibly  in  planning  replenishment  at  sea  operations, 
with  more  accurate  predictions  resulting  in  each  case. 

The  formulation  and  programming  of  a  simulation  model  based  on 
these  alongside  time  distributions,  and  thus  sensitive  to  changes  in 
load  requirements  is  the  second  objective.   While  incorporating  the 
above  known  features,  the  simulation  should  be  designed  to  enable  its 
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user  to  easily  adjust  the  arrangement  of  the  ships  involved,  load 
requirements,  and  distribution  parameters.   In  this  regard,  he  could 
utilize  it  to  test  the  suitability  of  various  arrangements  of  combatant 
ships  and  the  sensitivity  of  load  changes. 
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III.   SOURCE  OF  DATA 

Data  used  to  obtain  the  distribution  of  alongside  completion  times 
were  from  replenishment  at  sea  (RAS)  reports  3180/1A  and  3180/lB  for 
the  period  December  1967  to  April  1969.   Commander  Service  Force,  U.S. 
Pacific  Fleet  requires  these  reports  be  made  to  him  immediately 
following  underway  replenishment  by  those  replenishment  ships  under 
his  cognizance.   Each  report  includes  the  replenishment  ship  type 
and  hull  number,  receiving  ship  type  and  hull  number,  total  alongside 
time,  quantity  transferred  in  barrels  or  short  tones  as  appropriate, 
UNREP  speed,  whether  day  or  night,  whether  EASTPAC  or  WESTPAC,  and 
whether  or  not  there  was  simultaneous  port  and  starboard  UNREP  in 
progress.   Full  instructions  concerning  the  submission  of  these 
reports  are  available  in  COMSERVPAC  INST  3180. 3F  of  6  November  1967. 

Alongside  time  is  defined  as  the  elapsed  time  from  commencing 
approach  from  300  yards  astern  to  time  of  departure.   Only  data 
submitted  by  oiler  (AO)  and  ammunition  (AE)  type  ships  involving 
destroyer  (DD)  and  attack  carrier  (CVA)  recipients  were  used,  since 
a  sufficient  quantity  was  available  only  for  these  combinations. 
However,  sufficient  data  should  become  available  in  subsequent  years 
to  make  similar  analyses  possible  for  all  combinations  of  participants. 
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IV.   METHOD  OF  DATA  ANALYSIS 

In  order  to  accomplish  the  objective  of  a  more  accurate  estimation 
of  alongside  time  distributions  for  various  participants,  data  was 
first  categorized  as  to  type  of  ships  involved,  then  separated  with 
regard  to  day  or  night  operation,  and  finally  subdivided  into  sets 
based  on  quantity  transferred.   The  distinction  based  on  day  or  night 
conditions  was  considered  necessary,  since  from  personal  experience 
and  as  shown  by  Besecker,  night  operations  consistently  take  longer 
than  those  conducted  in  daylight  hours.   The  load  intervals  were 
selected  on  the  basis  of  quantity  of  data  available  and  uniformity 
of  interval  ranges.   The  latter  was  desirable  in  order  to  establish 
load  intervals  necessary  for  the  computer  simulation. 

Data  in  each  category  was  analyzed  by  constructing  histograms  of 
the  relative  frequency  with  which  the  alongside  replenishment  time 
fell  within  ten  minute  intervals.   The  last  interval  is  always  a  half 
open  interval  from  its  selected  lower  bound  to  infinity.   In  each  case 
the  histogram  suggested  that  the  distribution  of  the  completion  times 
for  the  respective  load  intervals  could  be  fitted  to  a  gamma  or 
exponential  distribution  presented  in  Appendix  A.   This  was  expected 
since  the  distribution  of  alongside  times  over  all  intervals  had  been 
fitted  to  a  gamma  distribution  in  Ll ]• 

Figure  (1)  is  the  histogram  constructed  for  the  case  in  which  an 
oiler  (AO)  was  transferring  1000-1099  barrels  of  fuel  to  a  destroyer 
(DD)  during  daylight  hours.   This  is  presented  as  a  typical  example. 
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It  was  appropriate  to  apply  a  Chi-Square  goodness-of-f it  test  to 

the  null  hypothesis 

H  :  P,  =  ?4      i  -  (l,...,k) 
o   l    io 

against  all  alternatives,  where  P.  is  the  probability  of  finding  a 
completion  in  the  i   interval  with  respect  to  the  distribution  in 
question,  P.   is  the  relative  frequency  of  the  data  falling  within  the 
same  interval,  and  k  is  the  number  of  intervals.   The  probability  P 
was  calculated  using  Besecker's  computer  subroutine  GAMDIS  with  the 
gamma  parameters  estimated  from  the  sample  mean  and  variance  using 
the  relationships 


-2 

A    - 

T 

r 

2 

s 

^ 

_ 

\- 

T 

2 
s 

2 
where  T  is  the  sample  mean  and  s   the  sample  variance. 

The  Chi-Square  test  consisted  of  calculating  the  statistic 

K 

n        V   (X.-nP.  )2 


^io 


i=  1 


where  X   is  the  number  of  datum  which  fell  within  the  i    ten  minute 
interval  and  n  is  the  sample  size.   In  accordance  with  statistical 
theory   3  ,  if  H   is  true,  then  the  random  variable  Q  has  an  approxi- 
mate Chi-Square  distribution.   The  hypothesis  was  accepted  if 

Q  <X  (-95,  k-3). 

The  degrees  of  freedom  was  necessarily  k-3  since  two  parameters  were 
estimated  from  the  data.   Many  authors  caution  that  nP  ,  i=l,2  ...   k 
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should  be  greater  than  five,  since  Q  has  only  an  approximate  Chi-Square 
distribution.   This  warning  was  strickly  adhered  to,  which  necessitated 
combining  two  or  more  time  intervals  when  the  quantity  of  data  in  a  load 
interval  was  small. 

In  the  cases  involving  attack  carriers,  the  hypothesis  was  accepted 
for  all  load  intervals.   Since  the  Erlang  distribution  is  more  ameanable 
to  analytical  and  computer  simulation  models  and  only  differs  from  the 
gamma  distribution  in  that  r  is  restricted  to  be  a  positive  integer,  each 
r  was  rounded  to  the  nearest  positive  integer,  the  parameter  "X  adjusted 
according  to  the  expression 

"X  -JL 

T 

and  a  similar  test  conducted.   Again,  the  hypothesis  was  accepted  in 
each  case.   Results  of  these  latter  tests  are  given  in  Tables  I,  II, 
III,  and  IV  in  Appendix  B.   Figure  (2)  shows  a  typical  example  with  the 
probability  density  function  superimposed  on  the  relative  frequency 
histogram  for  a  case  involving  AO-CVA  participants.   The  figure  illus- 
trates the  method  of  combining  intervals  consistent  with  the  greater 
than  five  rule  given  above. 

When  the  above  method  was  used  for  cases  involving  destroyer  type 
ships,  the  hypothesis  was  rejected  in  each  case.   Re-examination  of  the 
histograms  for  these  cases  suggested  that  a  shift  of  the  time  axis  by 
an  amount  a  would  lend  to  an  acceptance  of  the  gamma  distribution. 
Several  approaches  were  taken  to  estimate  the  distribution  parameters 
r,   ,  and  a.   Although  many  were  not  successful  in  this  study,  all  are 
presented  for  reference  in  future  data  analyses. 

Noting  that  the  first  three  centeral  moments  of  the  gamma  distribu- 
tion given  by 
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E  \_t\  =  -|+  a,  (1) 

E    (T  -/ic)    =  I-  , 

V 

and 

E    (T  -  O3   =2L 

A3 

could  be  used  to  estimate  the  parameters,  the  following  relations  were 

derived: 

a      6 

r  =  4^  (2) 

m3  ■ 

*       2 

)  =  2l  (3) 

m3 

2  =  T  -  2  ^  (4) 

m 
3 

where  T  is  the  sample  mean  given  by 


n 


T  -  1 


=  i  y  t. 


i=l 


2 
s   is  the  sample  variance  given  by 

n 


s2  =  -L-     \j     (T  -  T)2  , 
n~l      1-1    X 


and  m   ^s  the  third  centeral  moment  given  by 

n 


1  ">        "  3 
m.  =  -  L    (T  -  T)  . 

3   n   1=1  i 


By  rounding  the  shape  parameter  r  to  the  nearest  positive  integer  and 
adjusting  the  scale  and  shift  parameters,  the  following  expressions 
were  obtained: 
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6 

(I       S    x* 

in 
3 


_2  » 


JT„ 


and 


a   —      V7T 
a  =  T  -  s  \/  r 


where  (x)*  is  the  nearest  integer  to  x.   Using  this  method  in  several 

A 

cases,  r  was  found  to  be  unity  in  each  case  suggesting  an  exponential 
distribution.   Figure  (3)  shows  a  typical  case  with  the  exponential 
distribution  superimposed  on  the  histogram.   The  hypothesis  was  rejected 
for  these  cases  as  clearly  suggested  by  the  figure. 

Concern  over  the  poor  results  and  the  fact  that  observed  times  were 
less  than  a,  prompted  another  investigation  into  the  method  of  obtaining 
the  shift  parameter.   Hypothesizing  an  exponential  distribution  and 
constructing  the  likelihood  function 


L(A,a)  =  |  A  e  for  at  min  T 

\t\in\    i 


the  maximum  likelihood  estimator  for  a  was  obtained  as 


a  =  min  T 


This  Implied  that 


A- 


T-S 
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from  Eq.  (1).   However,  for  all  cases  examined  using  this  method,  the 
hypothesis  was  rejected  on  the  basis  of  the  Chi-Square  goodness-of-fit 
test  as  illustrated  in  Figure  (4). 

Concluding  that  rounding  to  obtain  an  Erlang  distribution  was 
inappropriate  when  0  -  ri  1,  an  attempt  was  made  to  use  the  maximum 
likelihood  estimate  for  a  under  the  assumption  of  the  gamma  distribution, 
The  likelihood  function  in  this  case  is 


n 


L(a)  = 


Tn  n        r-1  ">^  (Ti  "  a) 

~\       T\  (T  -  a)    e   1  ~  1      for  aimin  T* 


1  0  otherwise, 

resulting  in  the  maximum  likelihood  estimate  of  a  as 


A 

a  =  min  T   , 


again.   Then  referring  to  the  expressions  obtained  using  the  method  of 
moments  and  adjusting  the  other  parameters  to  a,  the  following  equations 
were  derived: 

r  =  (T-a)2 
s^ 

A      _ 

>  =  ?-a 

2 
s 

Gamma  distributions  suggested  by  this  third  method  did  not  fit  the  data 
as  shown  in  Figure  (5) . 

In  another  attempt  to  use  the  method  of  moments  to  estimate  the 
distribution  parameters,  no  adjustment  was  made  to  the  estimates 
given  by  Eqs.  (2),  (3),  and  (4).   But  this  also  proved  unsatisfactory 
as  illustrated  in  Figure  (6). 

Convinced  that  the  gamma  distribution  was  appropriate  from  the 
hnpe  of  the  histograms  based  on  ten  minute  intervals  from  zero, 
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extensive  probability  tables  were  generated  using  Besecker's  subordinate 
GAMDIS.   Then  by  scanning  these  tables  for  similar  patterns  of  the 
relative  frequency  ratios  of  time  intervals  compared  to  the  probability 
that  a  completion  time  would  be  within  the  same  interval,  an  Erlang 
distribution  was  selected  that  appeared  to  fit  the  data.   A  Chi-Square 
goodness-of-fit  te?t  was  conducted  and  the  hypothesis  H  accepted  if 


o 


2 
Q  X     (,95,k-l). 


This  method  proved  quite  fruitful,  resulting  in  the  acceptance  of 
the  hypothesis  associated  with  an  Erlang  distribution  in  each  case.   A 
typical  case  is  shown  in  Figure  (7) .   Results  of  these  tests  are  given 
in  Tables  V,  VI,  VII,  and  VIII  in  Appendix  B  with  T  calculated  using 
the  expression 

A        A        A 

T  =  JL  +  a. 

The  estimate  a  was  chosen  as  the  lower  bound  of  the  first  interval  in 
which  data  was  observed.   Figures  (12)  through  (56)  in  Appendix  C 
illustrate  the  fit  of  the  probability  density  functions  and  the  histo- 
grams for  each  case  involving  destroyer  type  ships. 
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V.   THE  COMPUTER  SIMULATION  MODEL 

As  stated  in  Nay lor  \_5\   the  purpose  of  simulating  a  large  operation 
such  as  replenishment  at  sea  is  to  conduct  situational  experiments  that 
would  ordinarily  be  too  expensive  and/or  cumbersome  to  perform  physi- 
cally.  Because  of  the  many  conceivable  combinations  of  participant 
ships  and  variations  in  quantities  transferred  between  them  that  one 
might  want  to  examine,  physical  experimentation  would  be  unreasonable 
here.   Therefore,  based  on  the  distributions  of  the  previous  section, 
a  Monte  Carlo  type  digital  computer  simulation  of  the  replenishment  at 
sea  operation  named  SIMA  was  devised  and  programmed  in  accordance  with 
the  second  objective. 

SIMA  was  designed  to  accommodate  a  system  of  from  one  to  four 
replenishment  ships,  composing  an  Underway  Replenishment  Group  (URG) , 
with  each  ship  having  one  or  two  replenishment  stations  on  the  port  or 
port  and  starboard  sides.   A  maximum  of  twenty  combatant  ships  can  be 
in  the  queue  at  each  station.   Furthermore,  six  types  of  combatant 
ships  can  be  in  the  system.   This  capacity  is  well  beyond  what  one 
would  normally  expect  during  an  actual  replenishment  operation, 
therefore  overload  of  the  model  is  not  expected. 

Drawing  from  personal  experience,  the  model  was  formulated  as 
realistically  as  possible,  thereby  taking  advantage  of  the  known 
factors  of  a  replenishment  at  sea  system.   All  combatant  ships  are 
initially  distributed  among  the  various  stations  of  the  supply  ships, 
Lth  each  individually  characterized  by  the  load  it  requires  from  each 
•  1  <ni  ■  hn.-nt  -hip.   Combatant  ships  proceed  sequentially  through  the 
ind  upon  arrivnl  at  a  ropl onishment  ship  it  enters  the  smaller 
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(if  more  than  one  exists)  queue.   Each  ship  receives  the  services  of 
every  supply  ship  and,  having  done  so,  is  ejected  from  the  system. 
Although  the  combatant  ships  do  pass  through  the  system  in  a  specified 
manner  from  their  initial  positions,  these  conditions  closely  resemble 
the  ability  of  a  Task  Group  Commander  to  formulate  the  operation  prior 
to  its  commencement,  utilizing  his  knowledge  of  the  approximate  loads 
transferred  among  the  participants. 

Using  the  simulation  as  an  experimental  device,  one  is  able  to 
specify  the  participants,  number  of  stations,  initial  sizes  of  the 
station  queues,  and  quantities  transferred  between  the  respective 
ships,  and  then  periodically  observe  the  state  of  the  system.   These 
observations  are  made  at  variable  time  increments  which  are  the  random 
alongside  completion  times  generated  as  per  Appendix  D.   An  observation 
is  made  and  the  system  advances  whenever  any  one  of  the  URG  ships  has 
completed  a  replenishment.   Each  observation  consists  of  the  time  since 
the  start  of  the  operation  and  the  state  of  the  system.   The  state  of 
the  system  at  these  times  is  the  number  of  ships,  including  the  one 
being  served,  at  each  replenishment  station  and  is  represented  by  an 
n  dimensional  vector,  where  n  is  the  number  of  stations  in  the  system. 

The  i   component  of  the  state  vector  being  the  number  of  ships  at  the 

th 
i   station. 

This  model  was  programmed  in  Fortran  IV  and  tested  on  the  IBM 

model  360-67  digital  computer  at  the  U.S.  Naval  Postgraduate  School, 

Monterey,  California.   A  printout  of  the  program  is  presented  In 

Appendix  E.   Because  of  the  digital  nature  of  the  simulation,  ship 

types  were  given  the  following  identification  numbers: 


U 


Oiler  -  1 

Ammunition  ship  -  2 
Destroyer  -  1 
Attack  carrier  -  2 

Since  all  the  parameters  derived  in  the  previous  section  are  given  in 

data  statements,  they  may  be  revised  by  changing  the  statements.   In 

using  this  technique,  a  user  need  only  supply  the  following  information 

to  conduct  replenishment  at  sea  experiments: 

1.  Type  and  sequence  of  supply  ships 

2.  Day  or  night  operation 

3.  Type,  initial  position,  and  load  requirement  of  the  combatant 
ships . 

This  information  is  provided  by  a  data  deck.   The  first  m(m  -  1,2,3,  or 

4)  cards  represent  the  m  supply  ships,  each  card  having  the  following 

format: 

Columns  1-3    Type  of  replenishment  ship,  i.e.  AO,  AE 
Columns  4-5    Number  of  stations  available.   (1  or  2) 
Columns  6-7    Supply  ship  identification  number 
Columns  8-10   Blank 
Column  11      Must  be  either  0,1,  or  blank. 

0  or  blank  implies  that  there  is  at  least 
one  other  replenishment  ship  in  the  system  and 
another  card  must  be  read. 

1  implies  no  more  supply  ships  are  in  the 
system  and  the  program  continues  to  the  next 
statement. 

Columns  12-80   Blank 

The  order  in  which  these  cards  appear  designates  the  sequence  in  which 

combatant  ships  pass  through  the  system.   The  m+lst  card  specifies  day 

or  night  operation  and  has  the  following  format: 

Column  1       Blank 

Column  2       Must  be  either  1  or  2. 

1  designates  day 

2  designates  night 
Columns  3-80    Blank 

The  remaining  cards  give  information  concerning  the  combatant  ships, 
with  the  cards  arranged  as  these  ships  initially  appear  in  the  system. 
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The  remaining  cards  give  information  concerning  the  combatant  ships, 
with  the  cards  arranged  as  these  ships  initially  appear  in  the  system. 
The  initial  condition  of  each  station's  queue  for  each  supply  ship  is 
given  in  accordance  with  the  order  of  the  supply  ship  cards,  with  the 
port  station  designated  first.   Each  card  represents  a  combatant  ship 
with  information  given  in  the  following  format: 

Columns  1-3      Ship  type,  i.e.  DD,  CVA 

Columns  4-5      Ship  identification  number. 

Columns  6-11     Load  required  from  first  supply  ship. 

Columns  12-17    Load  required  from  second  supply  ship. 

Columns  18-23    Load  required  from  third  supply  ship. 

Columns  24-29    Load  required  from  fourth  supply  ship. 

Column  30        Must  be  either  0,  1,  or  blank. 

0  or  blank  indicates  that  there  is 
at  least  one  other  combatant  ship  in 
the  queue  for  the  station  in  question 
and  another  card  must  be  read. 

1  indicates  that  no  more  ships  are 
initially  in  the  queue  in  question. 

Columns  31-80    Blank 

Figure  (8)  shows  a  hypothetical  initial  arrangement  of  the  supply 
and  combatant  ships,  illustrating  the  possible  routes  of  the  combatant 
ships.   As  an  example,  this  system  was  simulated  using  SIMA  in  day  and 
night  operation.   Figure  (9)  shows  the  data  cards  submitted  to  set  up  the 
system  for  the  daylight  case.    Since  only  two  supply  ships  are  in  the 
example,  the  results  may  be  represented  graphically  as  in  Figures  (10) 
and  (11).   The  origin  represents  the  starting  time  of  the  operation. 
Each  dot  represents  a  combatant  ship  having  completed  replenishment  at 
either  the  AO  or  AE.   The  length  of  the  line  segment  between  two 
consecutive  dots  indicates  the  elapsed  time  in  minutes  between  two 
consecutive  completions  of  the  replenishment  operation.   For  example, 
in  Figure  (10)  a  ship  departed  the  AO  at  approximately  60  minutes,  a 
ship  left  the  AE  at  approximately  80  minutes,  a  ship  left  the  AO  at 
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Figure  8.   Hypothetical  initial  arrangement  of  supply  and 
combatant  ships. 
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Figure  9.   Data  deck  used  to  set  up  example  problem. 


3S 


400 


300 


200 


100 


' 


♦  * 


0  —  i 

0      100 


200       300       400 
Ammunition  Ship  Time  (minutes) 


i 
500 


600 


Figure  10.   Graphical  representation  of  simulation 
results  for  the  example  system  during 
daylight  operation. 
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Figure  11 


Graphical  representation  of  simulation  results  for  the 
example  system  during  night  operation. 
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approximately  100  minutes,  a  ship  left  the  AE  at  approximately  100 
minutes,  and  so  on,  regarding  the  commencement  of  the  replenishment 
operation  0  hours  and  0  minutes.   Knowing  the  initial  arrangement  of 
the  ships,  one  is  able  to  reconstruct  the  movement  of  each  combatant 
ship  through  the  system  and  its  departure  time  from  the  supply  ships. 

Upon  running  the  simulation  fifty  times  for  the  example,  it  was 
found  that  the  expected  (average)  time  to  complete  the  replenishment 
operation  during  daylight  hours  was  598.62  minutes  while  night  operations 
required  635.38  minutes.   Figures  (12)  and  (13)  are  the  histograms  for 
the  frequency  ratios  of  total  replenishment  time  of  all  combatants  for 
these  experiments. 

We  see  from  the  data  cards  that  the  destroyer,  initially  the  second 
ship  in  the  queue  at  the  port  side  of  the  AO,  requires  648  barrels  of 
fuel  oil  and  2200  short  tons  of  ammunition.  When  this  was  changed  to 
1100  barrels  and  3000  short  tons  and  the  simulation  rerun,  the  expected 
time  to  complete  the  operation  was  605.54  for  day  and  638.10  for  night 
conditions.   Sensitivity  analyses  along  these  lines  could  easily  be 
accommodated  by  the  simulation  model. 

SIMA  might  also  be  useful  in  planning  replenishment  at  sea 
operations.   By  tracing  the  path  of  particular  ships  in  simulation, 
one  could  detect  possible  bottlenecks  in  a  proposed  replenishment 
system.   Then  by  repositioning  one  or  more  ships  in  the  formation  and 
conducting  additional  experiments,  one  could  plan  a  more  efficient 
arrangement  to  be  used  in  actual  replenishment  operations. 
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Relative  frequency  histogram  for  total 
replenishment  time  of  all  combatants  from 
fifty  simulation  runs  for  day  operation. 
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Figure  13.   Relative  frequency  histogram  for  total 

replenishment  time  of  all  combatants  from 
fifty  simulation  runs  for  night  operation 
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VI.   SUMMARY 

In  this  paper  we  have  shown  that  an  Erlang  or  exponential  distribu- 
tion may  be  fitted  to  the  CONREP  alongside  time  data  involving  AO,  AE, 
DD,  and  CVA  type  ships  for  each  of  several  load  intervals.   It  is 
possible  for  these  distributions  to  be  used  in  analytical  models 
employing  Laplace  transforms  that  would  be  sensitive  to  changes  in 
load  requirements  of  the  combatant  ships.   Furthermore,  several  methods 
for  obtaining  the  parameters  for  the  gamma  and  Erlang  distributions  are 
presented  for  reference  in  future  data  analyses. 

The  second  part  of  the  study  utilizes  the  distributions  in  a  compu- 
ter simulation  model.   Using  the  simulation  as  an  experimental  device, 
one  can  easily  and  realistically  investigate  the  suitability  of  various 
CONREP  situations  confident  that  the  results  will  reflect  one's  know- 
ledge of  the  approximate  loads  required  by  each  of  the  combatant  ships. 
It  was  demonstrated  that  SIMA  could  be  used  to  estimate  the  time  to 
complete  the  operation  and  investigate  the  sensitivity  of  certain  load 
requirements,  and  suggested  that  the  model  be  used  in  planning  underway 
replenishment  formations. 
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APPENDIX  A 
The  Gamma,  Erlang,  and  Exponential  Distributions 

The  gamma  probability  density  function  is  given  by 

f(x)  =     \   aV'V  x 

P(n)  for  xZ  0 

0  otherwise, 

where  r  is  the  shape  parameter  and  "X  the  scale  parameter.  In  shifting 
the  function  by  an  amount  a  the  random  variable  x  would  be  replaced  by 
x-a  in  the  above  expression. 

The  Erlang  and  exponential  distributions  are  special  cases  of  the 
gamma  distribution.   If  r  is  a  positive  integer,  the  probability 
density  function  may  be  written  as 

f(x)  =    \   ^   x   e 

(n-1)!  for  xZ   0 

0  otherwise, 

and  is  that  of  the  Erlang  distribution.   If  r  is  unity,  the  distribu- 
tion is  the  well  known  exponential  distribution. 
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APPENDIX  B 
Tables 

Table  I.   AO  vs .  CVA  -  day 


load  int. 

f 

r 

Q 

r< 

.95, k 

0-3999 

85.0 

.096 

8 

2.615 

5.99 

4000-4999 

98.7 

.244 

24 

2.716 

3.84 

5000-5999 

115.9 

.117 

14 

3.444 

5.99 

6000-6999 

127.2 

.074 

9 

1.497 

3.84 

7000-7999 

142.2 

.050 

7 

4.651 

5.99 

8000-8999 

139.3 

.156 

22 

2.411 

3.84 

9000-9999 

126.3 

.047 

6 

1.218 

5.99 

10000- 

143.3 

.060 

4 

2.584 

5.99 

Table  II.   AO  vs.  CVA  -  night 


load  int. 

T 

r 

Q 

X2( 

.95,k-3) 

0-5999 

118.2 

.068 

8 

1.007 

5.99 

6000-7999 

134.1 

.088 

12 

1.154 

5.99 

8000-9999 

139.9 

.063 

9 

0.856 

5.99 

10000-11999 

141.8 

.136 

19 

0.732 

5.99 

00- 

194.2 

.056 

11 

0.635 

5.99 
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Table  III.   AE  vs.  CVA  -  day 


load  int. 

T 

A 

A 

r 

Q 

t2 

C.95,k-3) 

0-9999 

74.7 

.092 

7 

2.486 

3.84 

10000-19999 

105.2 

.114 

12 

2.624 

3.84 

20000-29999 

159.9 

.111 

18 

3.055 

5.99 

30000-39999 

189.8 

.141 

27 

2.180 

3.84 

40000- 

229.4 

.079 

18 

0.067 

3.84 

Table  IV.   AE  vs.  CVA  -  night 


load  int. 

T 

A 

r 

Q 

2 

X 

(.95,k-3) 

0-19999 

90.1 

.084 

8 

2.573 

3.84 

20000-29999 

162.8 

.033 

6 

0.541 

3.84 

30000-39999 

168.8 

.076 

13 

1.402 

3.84 

40000-49999 

216.7 

.113 

24 

0.649 

3.84 

50000- 

258.1 

.094 

24 

0.823 

5.99 
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Table   V.      AO  vs.    DD  -    day 


A 

* 

A 

A 

X2(.95,k-1) 

load  int. 

T 

T 

"A 

r 

a 

Q 

0-399 

41.3 

40.0 

.100 

2 

20 

4.558 

11.10 

4Q0-499 

44.9 

50.0 

.100 

3 

20 

2.385 

9.49 

500-599 

49.8 

50.0 

.098 

2 

30 

1.973 

9.49 

600-699 

43.1 

40.0 

.100 

2 

20 

1.104 

9.49 

700-799 

47.7 

51.4 

.140 

3 

30 

3.172 

7.81 

800-899 

52.9 

53.3 

.060 

2 

20 

7.912 

12.60 

900-999 

52.8 

55.0 

.120 

3 

30 

7.143 

9.49 

1000-1099 

49.6 

55.0 

.120 

3 

30 

4.413 

14.10 

1100-1199 

53.0 

50.0 

.100 

3 

20 

5.529 

12.60 

1200-1299 

53.4 

50.0 

.100 

3 

20 

8.167 

12.60 

1300-1399 

56.4 

63.3 

.140 

5 

30 

1.389 

7.81 

1400-1499 

54.1 

55.0 

.040 

1 

30 

7.137 

12.60 

1500-1599 

62.6 

58.5 

.070 

2 

30 

4.637 

12.60 

1600-1699 

65.3 

65.7 

.140 

5 

30 

4.189 

12.60 

1700-1799 

61.3 

60.0 

.100 

3 

30 

3.306 

9.94 

1800-1899 

64.1 

60.0 

.100 

3 

30 

1.509 

11.10 

1900-1999 

67.3 

70.0 

.100 

3 

40 

2.439 

9.94 

0-2099 

67.4 

40.0 

.100 

3 

40 

3.073 

11.10 

2100-2199 

65.2 

65.0 

.040 

1 

40 

2.306 

9.49 

2200- 

68.3 

70.0 

.100 

3 

40 

4.925 

11.10 
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Table   VI.      AO   vs.    DD  -   night 


load  int. 

T 

A 

T 

A 

A 

r 

i 
a 

Q 

,'.C.95,k-l) 

0r599 

65.6 

67.5 

.080 

3 

30 

5.626 

11.10 

600-699 

63.3 

67.5 

.080 

3 

30 

1.510 

9.49 

700-799 

70.2 

67.5 

.080 

3 

30 

0.253 

7.81 

800-899 

78.6 

67.5 

.080 

3 

30 

1.784 

7.81 

900-999 

67.8 

67.5 

.080 

3 

30 

0.318 

7.81 

1000-1099 

61.3 

67.5 

.080 

3 

30 

1.313 

9.49 

1100-1199 

62.1 

67.5 

.080 

3 

30 

2.912 

7.81 

1200-1299 

56.1 

65.0 

.080 

2 

40 

1.799 

11.10 

1300-1399 

67.6 

65.0 

.080 

2 

40 

8.014 

9.49 

1400-1499 

71.6 

67.5 

.080 

3 

30 

1.528 

7.81 

1500-1599 

55.9 

65.0 

.160 

4 

40 

3.191 

7.81 

1600-1699 

74.4 

75.5 

.080 

3 

40 

0.549 

7.81 

1700-1799 

78.8 

75.5 

.080 

3 

40 

1.854 

7.81 

1800- 

76.7 

75.5 

.080 

3 

40 

4.563 

9.49 

AS 


Table  VII.   AE  vs.  DD  -  day 


load  int. 

T 

A 

T 

r 

A 

a 

Q 

X2 

(.95,k-l) 

0-999 

42.8 

46.0 

.250 

9 

10 

3.584 

9.49 

1000-1999 

63.3 

60.0 

.100 

4 

20 

1.013 

9.49 

2000-2999 

77.0 

80.0 

.080 

4 

30 

3.332 

12.60 

3000-3999 

97.8 

106.6 

.060 

4 

40 

7.927 

14.10 

4000-4999 

110.5 

105.4 

.044 

2 

60 

1.524 

11.10 

5000- 

127.5 

120.0 

.080 

4 

70 

5.642 

12.60 

Table  VIII.   AE  vs.  DD  -  night 


load  int. 

T 

A 

T 

A 

A. 

r 

A 

a 

Q 

X2( 

.95,k-l) 

0-999 

45.6 

40.0 

.100 

3 

10 

5.634 

9.49 

1000-1999 

72.3 

67.5 

.080 

3 

30 

1.075 

9.49 

2000-2999 

90.9 

87.5 

.080 

3 

50 

2.913 

9.49 

3000-3999 

104.7 

100.0 

.040 

2 

50 

3.666 

9.49 

4000- 

139.7 

145.0 

.080 

3 

70 

3.149 

12.60 
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APPENDIX  C 
Figures 

This  appendix  contains  figures  illustrating  the  fit  of  the 
probability  density  functions  to  the  histograms  for  all  cases 
involving  destroyer  type  ships. 
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APPENDIX  D 
Generating  Random  Observations  from  an  Erland  Distribution 

The  most  important  operation  in  a  Monte  Carlo  type  computer 
simulation  is  the  generation  of  appropriate  random  variates.   In  this 
study,  it  was  necessary  to  obtain  random  variates  from  the  Erland 
distribution.   Methematically  the  Erland  distribution  is  a  convolution 
of  r  independent  identical  exponential  distributions,  i.e.,  the 
distribution  of  the  sum  of  r  exponentially  distributed,  independent 
random  variables  whose  density  is  given  by 

f(y)  =  e~>y,  Y>0. 

Thus,  we  are  able  to  generate  Erlang  random  variables  by  summing  r 
independent  exponentially  distributed  variables  Y  ,  Y  ,  Y  ,  ...  ,  Yr, 

each  with  expected  value  1/   .   Doing  so  we  obtain  that 

r 
x  = 


I'i 


has  the  Erland  distribution  with  parameters  r  and   ,  i.e.  with  density 


_  r  r-1  -Ax 
>   x    e 
f(x)  = 


(r-1)  !  for  x  1  0 

0  otherwise. 

The  question  now  is  how  to  obtain  exponentially  distributed  random 

variates?   The  cumulative  distribution  function  for  the  exponent i«l 
distribution  is  given  by 
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f(y)  =  1-e""*7,  y7-0. 


It  is  well  known  that  F(Y),  and  by  symmetry  l-F(Y)  also,  is  a  uniform 
random  variable  over  the  interval  (0,1)  regardless  of  what  the  distri- 
bution of  the  continuous  random  variable  Y  is.   This  is  equivalent  to 
setting 

-*Y 
W  =  e 

where  W  is  a  random  variable  from  the  uniform  (0,1)  distribution. 
Therefore,  taking  the  natural  logarithm  of  both  sides 

-*Y  -  ln(W) 

which  yields 

Y  =  -  -L  ln(W). 
<\ 

Thus,  if  W  is  a  uniform  random  variable  over  (0,1),  Y  will  have  the 

desired  exponential  distribution.   Substituting  into  Eq.  (1)  we  get 


X  =-  -l-ln(TT  V- 


There  are  several  methods  available  for  generating  uniform  (0,1) 
random  variables.   The  IBM  Scientific  Subrouting  RANDU  was  used  in 
the  study. 
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APPENDIX  E 
Computer  Program 

This  appendix  contains  a  printout  of  the  computer  program  SIMA 
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